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Abstract 



We study a 3D generalization of the toric code model introduced recently by Chamon. This 
is an exactly solvable spin model with six-qubit nearest-neighbor interactions on an FCC lattice 
whose ground space exhibits topological quantum order. The elementary excitations of this 
model which we call monopoles can be geometrically described as the corners of rectangular- 
shaped membranes. We prove that the creation of an isolated monopole separated from other 
monopoles by a distance R requires an operator acting on U(R 2 ) qubits. Composite particles 
that consist of two monopoles (dipoles) and four monopoles (quadrupoles) can be described as 
end-points of strings. The peculiar feature of the model is that dipole-type strings are rigid, 
that is, such strings must be aligned with face-diagonals of the lattice. For periodic boundary 
conditions the ground space can encode 4g qubits where g is the greatest common divisor of 
the lattice dimensions. We describe a complete set of logical operators acting on the encoded 
f*"^ . qubits in terms of closed strings and closed membranes. 
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1 Introduction 



In the last decade, many two-dimensional spin models which exhibit topological order have been 
proposed and studied. Well-known examples include the 2D toric code and quantum doubles of 
finite groups p], Levin- Wen string-net models {2J, topological color codes [3], Kitaev's honeycomb 
model [I], quantum dimer models on the triangular lattice [3 |6], and quantum loop gases [7J. The 
corresponding topological phases can be described by a suitable class of anyons [I] which captures 
many essential features of a model including the superselection sectors, the ground state degeneracy, 
and the behavior of excitations under braiding and fusion. 

Much less is known about topological order in three dimensions and a general classification of 
topological phases in 3D largely remains an open problem. As a first step in this direction, a 3D 
generalization of 2D toric codes has been studied by several authors [U [9j [10] . The 3D toric code can 
be viewed as a 3D lattice Z2 gauge theory. In this theory electric charges are point-like excitations 
that can be created by string-like operators. On the other hand, magnetic excitations correspond to 
closed loops of Z2 flux. Such excitations can be created by membrane-like operators. The 3D toric 
code features a macroscopic energy barrier for the logical membrane-like X operator, but not for 
the string-like logical Z operator. This implies that one encode a classical bit in such system whose 
value will be protected at finite temperature [9]. 

Other 3D models with both string and membrane-like logical operators have been analyzed in [SJ 
ITTj . In these models membrane-like and string-like operators always appear in pairs, such that the 
two operators anticommute if the string crosses the membrane at an odd number of points. Hence, 
one may be led to believe that some kind of duality must be at play between the supports of the 
logical operators. Such duality would say, for example, that if one logical operator of a qubit, say X, 
is membrane-like, then the other logical operator, say, Z, must be string-like or even point-like. Such 
duality would first of all imply that the distance of 3D error-correcting codes would be bounded by 
the linear size of the system. Secondly, it would provide strong evidence for a no-go result concerning 
thermal stability of a 3D passive quantum memory [T2] . 

In this paper we analyze a spin model originally introduced by Chamon [13] in which some of these 
beliefs can be examined with a greater care. The elementary excitations of the model are point-like 
particles which we call monopoles. Monopoles carry a non-trivial topological charge meaning that no 
local operator can create an isolated monopole from the ground space. In contrast to all previously 
studied 3D topological spin models, monopoles cannot be created by string-like operators. Rather, 
monopoles can be described as the corners of rectangular-shaped membranes. More specifically, 
suppose an excited state contains an isolated monopole separated from other monopoles by a distance 
R. We prove that creating such a state starting from the ground space requires an operator acting 
on Q(R 2 ) qubits. It implies that local errors cannot induce diffusion in a dilute gas of monopoles 
since only highly non-local operators can move monopoles between adjacent sites [13J. Composite 
particles that consist of two monopoles (dipoles) and four monopoles (quadrupoles) can be described 
as end-points of strings. The peculiar feature of the model is that dipole-type strings are rigid, that 
is, these strings must be aligned with face-diagonals of the lattice. The quadrupole-type strings 
are partially flexible, meaning that such strings can only be deformed within planes orthogonal to 
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body-diagonals of the lattice. 

The degenerate ground space of the model defines a quantum error correcting code with k = 4g 
logical qubits, where g is the greatest common divisor of the lattice dimensions. We describe a 
complete set of logical operators in terms of closed strings and closed membranes for the special 
case g = 1. The rigidity of dipole-type strings in our model leads to very interesting features of the 
corresponding quantum code. Specifically, if the lattice has periodic boundary conditions and the 
lattice dimensions are pairwise co-prime, a rigid string aligned with a face-diagonal must fully fill up 
a two-dimensional plane before it gets closed. Hence closed rigid strings can become membrane-like 
objects. We propose a subsystem encoding of a single logical qubit in which the only relevant logical 
operators are those associated with closed rigid strings and closed membranes. In this encoding 
closed flexible (quadrupole-type) strings can only affect the gauge subsystem. Several strategies of 
minimizing the weight of logical operators composed of rigid strings are discussed. 

It is important to emphasize the essential difference between membranes describing logical op- 
erators in our code and the ones in the 4D toric code [H]. In our model membranes must have a 
rectangular shape and the membrane operator creates excitations only near the corners of the rect- 
angle. Such membranes can be increased in size without paying extra energy penalty. In contrast, 
membranes in the 4D toric code may take arbitrary shapes and the membrane operator creates ex- 
citations along the entire boundary of the membrane. Increasing the size of such membranes costs 
energy growing linearly with the length of the boundary. The presence of this "energy barrier" is 
responsible for the thermal stability of 4D toric code [T5] manifesting itself in the exponentially large 
(as a function of the lattice size) relaxation time of the error-corrected logical operators. Although our 
3D model does not feature a linearly growing energy barrier similar to the 4D toric code, its dynamics 
towards a thermal equilibrium might be characterized by glassiness in the creation/annihilation of 
isolated monopoles as was argued in [13]. Such glassiness might result from the fact that local errors 
cannot induce diffusion of isolated monopoles. We do not expect that the relaxation time of the 
error-corrected logical operators in our 3D model will grow with the lattice size, see the discussion 
in Section 14.31 

Our main technical results pertaining to the Chamon's model [13] which we introduce formally 
in the next section can be summarized as follows. 

• Complete classification of bulk excitations 

• Lower bound on the weight of operators creating isolated monopoles 

• General formula for the ground space degeneracy on a 3-torus 

• Subsystem encoding of a qubit using only closed rigid strings 

• Proof of the zero-temperature stability 
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Figure 1: The generator S u and two possible ways for a pair of generators to have overlapping 
supports. Black dots indicate qubit locations. Generators S u are centered at sites indicated by open 
dots. 

1.1 The spin model 

We consider a 3D cubic lattice A with periodic boundary conditions and linear dimensions L x ,L y ,L z , 
that is, 

A = Z Lx x Z Ly x Z Lz . 

We shall say that a site u = k) e A is even (odd) iff i + j + k is even (odd). Note that the parity 
of a site is well-defined for periodic boundary conditions only if all lattice dimensions L x ,L y , L z are 
even, since otherwise a transformation like % — >■ i — L x could change the parity of i + j + k. Let 
A even and A 0( m be the sublattices including all even and odd sites respectively. Note that each of the 
sublattices A even and A odd can be identified with face-centered cubic (FCC) lattice, see Fig. [2j We 
shall place a qubit at every even site of the lattice. Hence the total number of qubits is 



n 



— L x L y L z . 



We shall use the notations X u = a®, Y u = a^, Z u = o z u for the single-qubit Pauli operators acting on 
a site u G A even . Let x = (1, 0, 0), y = (0, 1, 0), and z = (0, 0, 1) be the basis vectors of the lattice. 
For every odd site u we define a stabilizer generator 

Su A" U _2; X U + X Y U — y Y U -\-y Z U ^ £ -^M+Z) (^) 

see Fig. [TJand Fig. [2J One can easily check that any pair of generators S u , S v commute, 

S U S V = S V S U for all u,v e A odd . 

Indeed, ciiiy p£iir S u ciiici S v may only overlap on a single qubit or on some pair of qubits, see Fig. [TJ 
In the first case S u and S v act on the shared qubit by the same Pauli operator and thus commute. 
In the second case S u and S v anti-commute at each of the two shared qubits, and thus overall they 
also commute. 

The generators S u are used to define a local term-wise commuting Hamiltonian [13] 



H 



u€A odd 



(2) 
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Figure 2: (Color Online) Qubits (black dots) live on the FCC lattice A even . The stabilizer generators 
(red) S u are centered on the open dots u in A odd . The double (blue) lines are examples of the six 
face- diagonals and the two dashed (green) lines are examples of the four body- diagonals. 

This model can be viewed as a natural 3D extension of the 2D toric code Hamiltonian. We illustrate 
the correspondence between the two models in Fig. [3j Our first result concerns the ground state 
energy and the ground space degeneracy of H. Note that the pairwise commutativity of the generators 
S u by itself does not guarantee that one can minimize all terms in H simultaneously, for example, if 
some product of the generators gives —I. Before stating the result let us introduce some notations. 
For any integers p, q, r we shall use notation gcd(p, q, r) for the greatest common divisor of p, q, r, 
that is, the largest integer that divides each of p,q,r. We shall parameterize the lattice dimensions 
as 

L x — 2p x , L y = 2p y , L z = 2p z . 

Define the subspace 

£ = m : S u |V) = |V) for all u E A odd }. 
In Appendix [A] we prove the following: 

Theorem 1. The subspace £ has dimension 2 Ag , where g = gcd(p x ,p y ,p z ) . 

Since states from £ minimize every term in the Hamiltonian, we conclude that £ is the ground 
subspace of H. Theorem [T] shows that the ground space degeneracy of the model Eq. fl2]) is determined 
not only by the topology of the underlying manifold (which is always the 3-torus), but also by the 
particular lattice discretization of this manifold. Hence this model is not purely topological in the 
same sense as the 2D toric code or the 3D models studied in [8]. It is worth pointing out that the 
choice of lattice discretization (even/odd lattice dimensions) does affect the ground state degeneracy 
for a slightly modified 2D toric code model studied by Wen [TB]. One reason why we call the model 
Eq. d2J) topological is that its ground state degeneracy is characterized by exponential insensitivity 
to weak local perturbations while the spectral gap above the ground state is stable against such 
perturbations. We can prove this by explicitly checking the sufficient conditions for stability of 
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(a) (b) (c) 

Figure 3: (a) The restriction of a generator S u onto a horizontal plane, (b) The correspondence 
between even sites lying in the chosen plane and links of the toric code lattice, (c) The change of 
basis X — y Y, Y — > X, Z — > —Z applied to every horizontal link maps the in-plane generators 
S u to the plaquette (X-type) and star (F-type) operators of the toric code. A string-like operator 
comprised only of Pauli Z's can commute with both plaquette and star operators only if it follows the 
diagonal of the lattice (the pink shaded region). Strings that contain X or Y Paulis fail to commute 
with the generators S u located in the adjacent horizontal planes. Hence extending the toric code to 
3D leads to rigidity of strings. Double circles indicate excitations created near the end-point of the 
string. 

topological phases derived in jTTj , see Appendix [B] In particular, we show that the ground subspace 
£ is a quantum error correcting code with a distance growing at least linearly with the lattice size. 

Our second result is a complete classification of bulk excitations of the model, see Sections [2] and [3) 
Since this classification does not depend on the boundary conditions, it is more natural to work with 
the infinite lattice, that is, A = Z x Z x Z. For the infinite lattice the ground state of the model is 
specified by the eigenvalue equations S u = 1 for all u G A 0( ^. The ground state is non-degenerate 
in the sense made precise in Appendix [Bl Since all generators S u pairwise commute and S 2 = I, 
one can describe any excited eigenstate of the model by specifying the eigenvalues S u = ±1 for each 
generator. The elementary excitations that we call monopoles flip the sign of a single generator. The 
following theorem proved in Section 13.11 shows that creating an isolated monopole (or any excited 
spot that contains an odd number of monopoles) separated from other excitations by distance R 
requires an operator acting on roughly R 2 qubits. 

Theorem 2. Let M C K odd be any finite subset containing odd number of sites. Let {Pr}r>i be a 
sequence of operators such that Pr creates excitations at every site u G M and may be some other 
excitations separated from M by distance at least R. Then Pr must act on Q(R 2 ) qubits in the limit 
R — y oo. 

The theorem immediately implies that the parity of the number of monopoles defines a Z 2 topo- 
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logical charge, that is, this quantity cannot be changed by any local operator. In this respect 
monopoles are similar to the electric and magnetic excitations in the 2D toric code. What is more 
surprising, Theorem [2] implies that no string-like operator can create an isolated monopole at its 
end-point. Indeed, otherwise one would be able to construct a sequence of string-like operators {Pr} 
as above whose weight grows only linearly with R. The absence of string operators associated with 
monopoles also implies that no local operator can create an isolated pair of monopoles from the 
ground state even if the two monopoles are separated by a constant distance^. Indeed, otherwise 
one would be able to create a long-range pair of monopoles separated by distance R by combining 
O(R) short-range pairs, which would require an operator acting only on O(R) qubits. It also shows 
that no local operator can induce hopping of monopoles between adjacent sites unless there are other 
monopoles nearby. In Section 13.11 we show that isolated monopoles can be created on the corners 
of rectangular-shaped membranes, see Fig. [TUJ It demonstrates that the lower bound in Theorem [2] 
is tight. The presence of isolated monopoles which can only created by membrane-like operators 
is an essential and new feature of the model. As was discussed in [13], it might have interesting 
physical implications, for example, the diffusive dynamics of the dilute gas of monopoles might be 
characterized by an exponentially small (as a function of the inverse temperature) diffusion constant. 

Let us now consider an excited spot M C K odd that contains an even number of monopoles. In 
Section 13.21 we show how to decompose M into a combination of composite excitations that we call 
dipoles and quadrupoles. These excitations are composed of two and four monopoles respectively. We 
show that both dipoles and quadrupoles can be created by string-like operators; in this respect they 
are similar to electric charges in the 3D surface code [8]. The unusual feature of the model is that 
dipoles can only be created at end-points of rigid strings, that is, these strings must be aligned with 
face- diagonals of the lattice. An example of a rigid string lying in a horizontal plane and aligned 
with the face-diagonal x + y is shown in Fig. [3]^c), see also Fig. [H A rigid string creates a dipole near 
each of its end-points. The quadrupole-type strings must lie in a pair of adjacent planes orthogonal 
to some body-diagonal of the lattice. We shall refer to such pair of planes as a bilayer. As we explain 
in Section [2TB] one can identify qubits of any bilayer with a 2D hexagonal lattice while the generators 
S u centered inside the bilayer can be identified with the hexagonal plaquette operators of Kitaev's 
honeycomb lattice model [4J. Taking a product of elementary link operators of Kitaev's model over an 
arbitrary path on the hexagonal lattice one obtains a string operator creating a quadrupole near each 
end-point of the chosen path. We call such strings flexible bilayer strings, or simply flexible strings, 
since they can be arbitrarily deformed as long as a string does not leave the bilayer it belongs to. In 
Section [3T21 we introduce a complete set of topological charges which we use to show that monopoles, 
dipoles, and quadrupoles are topologically distinct from each other. 

Properties of the quantum error correcting code corresponding to the model Eq. (j2j) are described 
in Section HJ Recall that the main parameters of a quantum code are the number of physical qubits 
n, the number of logical qubits k, and the code distance d which is the number of single-qubit errors 
needed to destroy the encoded information, see Section H] for formal definitions. Theorem [1] implies 

'"Note however, that the number of monopoles modulo four does not define a topological charge; for example, a 
single-qubit Pauli operator can transform 1 monopole into 3 monopoles |13j . 
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that the degenerate ground space C can be used to encode k = Ag logical qubits into n = Ap x p y p z 
physical qubits. The best lower bound on the distance d that we can rigorously prove is d — £l(L), 
where L is the smallest of the lattice dimensions, see Lemma |6] in Appendix [B] This lower bound 
however completely ignores some subtle features of the model such as the relationships between prime 
factors of the lattice dimensions. These features affect the ground state properties in a dramatic way 
as can be seen from Theorem [TJ Hence we believe that a more favorable scaling of the distance can 
be achieved by fine-tuning of the parameters p x ,P y ,Pz, see the discussion in Section [U 
We conclude by discussing some open problems in Section [5j 

2 Rigid and flexible strings 

In this section we describe two classes of string-like operators capable of creating remote clusters of 
excitations located near the two end-points of a string. Since classification of string operators is a 
bulk property we shall only consider an infinite lattice. 

The strings from the first class which we call rigid strings are straight lines aligned with one of 
the six face-diagonals of the lattice. The corresponding string operator creates a pair of excitations 
(a dipole) located near each end-point of the string. The strings from the second class which we call 
flexible strings lie in a pair of two adjacent planes orthogonal to one of the four body-diagonals of 
the lattice. We call such a pair of planes a bilayer. Flexible strings can follow an arbitrary trajectory 
within the chosen bilayer but they cannot leave it. The corresponding string operator creates four 
excitations (a quadrupole) at each end-point of the string. 

We shall see that string operators associated with closed flexible loops can be expressed in terms 
of the generators S u . Although a single rigid string cannot be closed into a loop, we shall see that one 
can combine multiple rigid strings into a closed 3-valent graph (a string-net) such as a tetrahedron. 
The corresponding string-net operator can also be expressed in terms of the generators S u . 

2.1 Notations 

Let us remind the reader of some standard notation pertaining to the FCC lattice. The basis vectors 
of the 3D simple cubic lattice A = Z x Z x Z will be denoted as 





' 1 " 




' " 




' " 


X = 





, V = 


1 


, z = 
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Let a,b,c G {1, 0, —1 = 1} be a triple of integers. We shall use the following abbreviations. 



[abc] 


a vector ax + by + cz 


[abc] -plane 


a plane orthogonal to the vector [abc] 
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In this notation the six face-diagonals of the lattice, see Fig. [2j are 



[110], [110], [101], [101], [Oil], [Oil], 



(3) 



and the four body-diagonals of the lattice are 



[111], [111], [111], [111]. 



(4) 



2.2 Rigid strings and dipoles 



Let m > be any integer. Consider a set of m + 1 sites 



7 = {(0,0,0), (1,1,0) 



. ., (m,m, 0)} C A 



even 



and an operator 



One can easily check that W(7) anti-commutes with exactly four generators S u for which u lies in 
the [001]-plane with the ^-coordinate k = and has exactly one neighbor in 7, see Fig. HJ namely, 
u = (—1, 0, 0), u = (0, —1, 0), u = (m + 1, m, 0), and u = (m, m + 1,0). The set of sites 7 is an 
example of a rigid string. 

The other rigid strings associated with the [010]-planes and [100]-planes can be obtained from 
this example by applying the lattice symmetries and replacing Z u by Y u and X u respectively. Let 
h be one of the six face-diagonals, see Eq. ([3]) or Fig. [2j We can formally define the rigid strings of 
type h as follows. 

Definition 1. A sequence of sites 7 = («o,Wi, • • • , u m ) with Ui G A even is called an rigid string of 
type h iff Ui + \ — Ui = h for all i — 0, . . . , m — 1. The sites uq and u m are called the end-points of the 
string. 

We shall refer to a pair of excitations associated with an end-point of an rigid string of type h as 
a dipole of type h. 

We may ask whether it is possible to define a closed rigid string 7 such that the corresponding 
string operator ^(7) is a product of the generators S u (similar to closed loops in the 2D toric code). 
Due to the rigidity of the strings one cannot close a string in a single plane. However one can form 
a three-dimensional object, a tetrahedron, by using all six types of rigid strings interconnected with 
each other to form a closed string-net, see Fig. [51 Note that the action of the three rigid strings at the 
qubits located at vertices of the tetrahedron cancels since XYZ = il. Thus the smallest tetrahedron 
operator in a cube of dimensions 3x3x3 equals the single generator S v at the center of the cube, 
depicted in Fig. [2j 

Here is how we define a general tetrahedron operator. Let C C A be a some cube of the lattice 
such that all eight vertices of C belong to A even . Let T be a tetrahedron formed by four vertices of 
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Figure 4: Example of a rigid string 7 lying in the [001]-plane. The string operator W{^f) acts by Z 
on the qubits in the shaded region. Double circles near the end-points of the string indicate locations 
of excitations created by the string operator Wi^y). A pair of excitations located near each end-point 
of the string is called a dipole. 
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Figure 5: (Color Online) Construction of a closed rigid string-net or tetrahedron operator. The six 
edges of the tetrahedron T are rigid strings of the six possible types lying on the faces of the cube. 
The string- net operator W(T) is the product of the six corresponding string operators. Note that 
W(T) acts trivially at the qubits located at vertices of T since the triple of strings incident to any 
vertex of T cancel each other through the identity XYZ ~ /. The open dot with two incident blue 
lines indicates the location of excitations created by a pair of string operators. The two excitations 
cancel each other. 
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Figure 6: (Color Online) Open dots indicate centers of generators S u located in the plane Si. Gray 
and black dots indicate qubits located in the planes E and E 2 respectively. The union of the planes 
E U E 2 is called a bilayer (shaded green area). A generator (red lines) centered in the plane Ex 
can be viewed as six-body hexagonal plaquette operators acting on this bilayer. The double circles 
indicate centers of generators S u located in the planes E_i and E 3 . These generators touch three 
qubits in E or three qubits in E 2 and can be represented as star operators, see Figs. [7] and [HJ 

C such that edges of T are diagonals of faces of C, see Fig. By construction, all sites lying on the 
six edges of T are even. We can denote the six edges of T as 7f , 7f> 7i, 7 2 > 7i, 7fj where 7" and 7 2 
is the pair of edges orthogonal to the axis a. We define the tetrahedron operator W (T) as 

w (?) = n n n <■ 

a=x,y,z j=l,2 

Each rigid string operator involved in W(T) creates a pair of excitations near the end-point of the 
string, i.e., near some vertex of the tetrahedron T. However, as one can easily check, the excitations 
created by a triple of strings incident to any vertex of T cancel each other and thus the tetrahedron 
operator W(T) commutes with all generators S v . A simple inspection shows that W(T) can be 
represented as a product of all generators S v in the interior of T (up to a phase factor): 

wcr)~ n ( 5 ) 

v£TnA odd 

The tetrahedron operators will be important in Section 13.11 Lemma [2J where we use them to detect 
an isolated monopole. 

2.3 Flexible strings and quadrupoles 

The construction of flexible strings is more involved. Let us start by constructing such a string lying 
in a pair of [lll]-planes So, E 2 , where 

E Q = {(i,j,k) e A : i + j + k = a}. (6) 
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Figure 7: The correspondence between the union of the sublattices EoUEiUE 2 and the 2D hexagonal 
lattice. The sites indicated by open dots (centers of hexagons) represent S x . The sites indicated by 
gray and black dots represent E and E 2 respectively. The generators S u located at Ei can be viewed 
as plaquette operators associated with the hexagons as shown on the right. 

We identify the union E U E 2 with the 2D hexagonal lattice f2h ex by projecting the sublattices E , 
E 2 onto some fixed [lll]-plane, see e. g. Fig. El We shall refer to the pair of planes E U E 2 as a 
[lll]-bilayer, or simply a bilayer. 

It will be convenient to color the sites of f2h ex in black and gray using the checkerboard coloring 
as shown in Fig. Then we have the following correspondence between the union of the sublattices 
Eq U Ei U E 2 and the 2D lattice Qhex'- 



3D 


2D 


E 


gray sites 


Si 


centers of hexagons (white sites) 


s 2 


black sites 



Define a triple of vectors 

\ x = [Oil], X y = [101], \ z = [110]. (7) 

Consider any pair of sites u G E and v G E 2 . By construction, (u, v) is a link of the hexagonal 
lattice iff v = u + A" for some a G {x, y, z}. We shall say in this case that (u, v) is an a-link. 
Define the link operator (as in Kitaev's honeycomb model [1]) 

K(u,v) = a Z a v if ( M 5 v ) i s an Qf-link. 

Given any pair of links e^e' the operators K e and K e i commute iff e, e' do not overlap and anti- 
commute iff e and e' share exactly one vertex. For any plaquette (hexagon) p and and any site w of 
fihcx define plaquette and star operators 

B p =\[ K e , A w = i J] K e , (8) 

e£dp eGstar(ui) 
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Figure 8: Plaquette and star operators. 



see Fig. [SJ Here dp is the set of six links forming the boundary of p, while star(iu) is the set of three 
links incident to w. Note that the action of A w onto the site w is proportional to iX w Y w Z w ~ /, 
that is, A w acts only on the nearest neighbors of w. In particular, A w acts only on gray (black) 
sites iff w is black (gray), see Fig. [HJ The only generators that can act non-trivially on qubits of 
i\cx = S U S 2 are those located in the planes S_l, Si, and S 3 . One can easily check that the action 
of these generators on the qubits of the hexagonal lattice can be described as follows: 



3D 


2D 


£-1 


star operators centered at black sites 


Si 


plaquette operators 


s 3 


star operators centered at gray sites 



We have the commutation rules 

[K e , B p ] = for all links e, for all plaquettes p. (9) 

Furthermore, 

[K e , A w ] = iff e is incident to w or e does not overlap with star(w). (10) 

Alternatively, 

K e A w = —A w K e iff e is incident to a neighbor of w but not incident to w. (11) 

Consider any pair of sites u, v G fihex (black or gray) and let 7 be any path on £\ex connecting u and 
v. Define a string operator 

win) = n K *- 

It follows from Eq. (Q that W^) commutes with all plaquette operators B p . In addition, since 7 has 
even number of links incident to any site w 7^ u,v, we infer from Eq. (fTTl) that W^) commutes with 
star operators A w unless w is located near the end-points u or v. We shall be mostly interested in 
the case when the end-points u, v are sufficiently far apart such that the commutation rules between 
Win) an d star operators can be analyzed independently near each end-point of 7. Then one can 
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Figure 9: Example of a flexible bilayer string. Double circles near each end-point of the string 
indicate the locations of excitations created by the string operator W{^). Recall that excitations 
located at gray and black sites correspond to generators S u with u G S 3 and u G £_i. The group of 
four excitations located near each end-point of the string is called a quadrupole. 

easily check that W (7)^ = — A W W(^) iff w is the end-point of 7, or w is a neighbor of the end-point 
of 7. 

It is clear that similar strings can be constructed in the bilayers orthogonal to any body-diagonal 
t, see Eq. (jlj). Such a flexible bilayer string lies in a pair of even planes orthogonal to some body- 
diagonal t separated by a single odd plane as Fig. [6] demonstrates. The group of four excitations 
located near each end-point of the string will be referred to as a quadrupole of type t. More formally, 
we define these string operators as 

Definition 2 (Flexible bilayer strings). Let t = (t x ,t y ,t z ) be any body- diagonal. A sequence of 
sites 7 = (w , Ui, . . . , u m ) with Ui G A even is called a flexible [t]-bilayer string iff there exists a sequence 
a , . . . , a m _i G {x, y, z} and e = ±1 such that Uj + i = Uj + (— l) J 'e X(aj) for all j = 0, . . . , m — 1, 
where the vectors A (a) are defined as 
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The sites Uq and u rn are called the end-points of the string. We can define a string operator W^) 
associated with 7 as 

m—l 

W(j) = H<T%a% +1 . (12) 

j=0 

We assume that the separation between the end-points of 7 is sufficiently large so one can derive 
the commutation rules near each end-point independently. The commutation rules between W(7) 
and the generators S v are then summarized by the following Lemma [TJ 
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Lemma 1. Let 7 = (1*0,1*1, . . . ,u m ) be a flexible [t]-bilayer string. Then W{~f)S v 
is one of the eight sites listed below: 
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Again, we can ask how one can construct closed flexible strings. It is not hard to see that such 
closed strings which lie, say, in a [lll]-bilayer, can be constructed from the hexagonal plaquette 
operators B p . Specifically, let p be any plaquette of flhex and v G Si be the center of p. By definition 
of the plaquette operators we have B p = S v and B p can be regarded as a string operator, B p = W(j) 
where 7 = dp is the closed flexible string that consists of the six links lying on the boundary of p. 
Similarly, if 7 is a closed flexible string without self-intersections, and Int(7) is the set of plaquettes 
encircled by 7, we get the identity 

W(7)~ J] B p- ( 13 ) 

pGlnt(7) 

We will use these closed flexible strings in Section 13.21 to characterize the topological charges. 



3 Classification of excitations 

The purpose of this section is to classify the excitations of our 3D model. We will first consider how 
single isolated excitations which we call monopoles can emerge. We say that an excited state 
has an isolated monopole at a site u 6 A odd iff S u \ip) = and S v \ip) = \ip) for all sites «/« 

in a sufficiently large neighborhood of u. Note that the string operators described in Section |2] can 
only create an even number of excitations (a dipole or a quadrupole) near each end-point of a string. 
Therefore a natural question is whether some more complicated string-like operators could create a 
single monopole. In Section [37X1 we answer this question in the negative: we prove in Lemma [2] that 
the minimum weight of a Pauli operator creating an isolated monopole grows as Q(R 2 ), where R is 
the separation between the monopole and the nearest excitation. We then use the lemma to prove 
Theorem |2j 

The goal of Section 13.21 is to identify all possible excitations of the model that are topologically 
non-trivial, that is, excitations that cannot be created locally from the ground state. We describe 
a complete set of Z 2 topological charges such that an excitation is topologically non-trivial iff it is 
characterized by non-zero value of some topological charge. 

In contrast to 2D systems that can always be described by a finite number of topological charges, 
3D systems may possess an infinite number of topological charges. The simplest example of such 
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Figure 10: Double circles indicate locations of four isolated monopoles created by a membrane 
operator. The membrane M consists of all even (black) sites that belong to the shaded region. The 
corresponding membrane operator is W(M) = YlueM 



behavior is a 3D stack of non-interacting 2D toric codes, in which each copy of the toric code 
contributes a constant number of charges. We shall see that in our model the number of topological 
charges is also infinite, but it becomes finite after factoring out quadrupole excitations. 

A simple corollary of our classification is the proof that rigid strings and flexible bilayer strings 
are the only possible strings in the theory. More precisely, let P be a Pauli operator creating an 
excited spot M C A 0( ^ with m = 0(1) excitations separated from all other excitations by a distance 
R 3> 1 (we are interested in the limit R — » oo). If m is odd, Theorem [2] shows that \P\ = Q(R 2 ), that 
is, P cannot be reduced to a constant number of string-like operators for any reasonable definition of 
a string. On the other hand, if m is even, we will prove that M can be decomposed into a constant 
number of dipoles and quadrupoles, see Section 13.31 and thus M can be created using a constant 
number of rigid strings and flexible strings. 



3.1 Monopoles and membrane operators 

Isolated monopoles can be created by operators supported on rectangular-shaped membranes such 
that there is one monopole sitting near each of the corners of the membrane, see Fig.[10j For example, 
consider a set of sites 

M = {(i,j,0) E A even : \i\ + \j\ < R} (14) 

and let 

W(M) = 11 Z u . (15) 

Assume for simplicity that R is even. Then one can easily check that W(M) creates four isolated 
monopoles at the sites listed below. 

W(M)S V — —S V W{M) iff v e {(R + 1,0,0), (0,i?+l,0), (-R— 1,0, 0), (0,-^-1,0)}. 

(16) 

The operator W (M) is an example of a membrane operator of weight 0(R 2 ). We shall now prove 
that the weight of any operator creating an isolated monopole must grow quadratically with R. 
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To illustrate the idea of the proof, let us first consider a 2D analogue of this problem. Suppose a 
Pauli operator P creates an electric charge e in the 2D toric code separated from other excitations 
by a distance R. The charge e can be detected by moving a magnetic charge m over an arbitrary 
closed loop 7 encircling e. The corresponding string operator W^) must anti-commute with P and 
hence the supports of P and W^) must share at least one qubit. Choosing a family of m = Q(R) 
pairwise disjoint loops 71, . . . ,7 m of increasing radius that encircle e, we conclude that the support 
of P shares at least one qubit with each loop 7$, and thus \P\ > m — Q(R). In our 3D model we 
shall detect isolated monopoles using a family of Q(R 2 ) closed rigid strings or tetrahedron operators, 
see Section 12.21 Since a tetrahedron has a one- dimensional support, we can choose Q(R 2 ) pairwise 
disjoint tetrahedra of linear size R enclosing the monopole to be detected (in the proof of the Lemma 
the tetrahedra are in fact not disjoint but overlap on 0(1) qubits). It yields a lower bound Q(R 2 ) 
on the weight of the operator P creating the monopole. Note that using operators with membrane- 
like support to detect the monopole would not be good enough, since we could only construct only 
0(R) non-overlapping membranes which enclose the monopole. Similarly, using closed flexible loops 
to detect the monopole would not be good enough, since one can only construct 0(R) encircling, 
non-overlapping, flexible loops in the bilayer that contains the monopole. 

Lemma 2. Let P be a Pauli operator such that PS U = —S U P for some site u G A 0( ^. Suppose 
PS V = S V P for all sites v 7^ u within distance R from u. Then \P\ > cR 2 for some constant c. 

Proof. Suppose one can construct m = Q(R 2 ) operators 0\, . . . , O m with the following properties: 

• Each operator Oj anti-commutes with P 

• Each qubit is acted on by at most c = 0(1) operators Oj 

Note that each qubit in the support of P can be responsible for anti-commutativity with at most c 
operators Oj. Since P anti-commutes with all m operators Oj, we conclude that c\P\ > m, that is, 
\P\ = Q(R 2 ). Let us explain how to choose the desired operators Oj. Let B be a ball of radius R 
centered at the site u occupied by the monopole. For any tetrahedron operator W(T) constructed 
as in Section [2721 such that u G T C B, one has the desired anti-commutation rule: 

PW{T) = -W{T)P for all T C B such that u G T. (17) 

Here we used the fact that W(T) can be represented as the product of all generators S u in the interior 
of T. It remains to choose sufficiently many tetrahedra with the right properties. Let T be a set of 
tetrahedra T of size R/10 such that mgT (and thus T C B). Obviously, \T\ = Q(R 3 ). Let us say 
that a pair of tetrahedra T,f e T is bad iff some edge of T and some edge of T' lie on the same 
line (and hence they can possibly overlap on Q(R) qubits). Otherwise, we shall say that T,T' is a 
good pair. Note that the supports of a good pair of tetrahedra T and T' overlap on the intersection 
of non-parallel edges, hence on 0(1) qubits. For any fixed T G T there are at most 0(R) tetrahedra 
T'gT such that the pair T, T' is bad. It follows that there exists a subset of m = Q(R 2 ) tetrahedra 
Ti, . . . , T m G T such that any pair T i; Tj is good: we can construct such a subset by picking a first 
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tetrahedron T\ , eliminating all tetrahedra which are bad for T\ and picking the next tetrahedron from 
the remaining set etc. At every step one looses O(R) tetrahedra, so starting with Q(R 3 ) tetrahedra 
allows us to pick m = Q(R 2 ) pairwise good ones. Choosing Oj = W(Tj), j = 1, . . . , m concludes the 
proof. □ 

Theorem [2] is a straightforward generalization of the above lemma. Indeed, suppose first that Pr 
is a Pauli operator. If R is much larger than the size of the excited spot M, one can use tetrahedra 
of size R to construct operators W{T) satisfying Eq. (fPTj) and thus prove the bound \Pr\ — Q(R 2 ). 
Suppose now that P is an arbitrary (non-Pauli) operator. Let be the ground state and \<f>) = P \ip) 
be an excited state with an isolated monopole, that is, S u \<p) = —\4>) and S v \4>) = \(j>) for all sites 
v 7^ u within distance R from u. Decomposing \<p) in the common eigenbasis of the generators 
S v , v G A odd , one can find an eigenstate \ip') such that {ip'\P\ip) ^ and \ip') contains an isolated 
monopole at u. Expanding P in the Pauli basis we conclude that at least one Pauli operator P' in 
this expansion must satisfy {il)'\P'\ip) ^ 0. Then by Lemma [2] one must have \P'\ = Q(R 2 ) and thus 
P acts on Q(R 2 ) qubits. 



3.2 A complete classification 

An excited eigenstate of the Hamiltonian Eq. (TSJ) can be described by a syndrome s that assigns 
an eigenvalue (— to any generator S u . We shall only consider excited states with finite energy 
(recall that in this section we consider an infinite lattice). More formally, let V be the group of 
Pauli operators (with possible infinite support). A function s: A odd — > {0, 1} is called a syndrome iff 
s(u) = 1 only for finitely many u and there exists a Pauli operator P G V such that 

S U P = (-1) S ^PS U for all u G A odd . (18) 

We say that P causes the syndrome s. Note that the commutation between S u and P is well-defined 
even if P has infinite support (although the commutation between Pauli operators with infinite 
support is not well-defined). One can also regard syndromes as binary strings with a finite Hamming 
weight whose bits are labeled by sites u G A odd . A bitwise sum of strings s, s 1 will be denoted s © s'. 
We shall group syndromes into equivalence classes called superselection sectors using the following 
definition. 

Definition 3. Two syndromes s and s' are equivalent iff s® s' is caused by a Pauli operator P G "P 
with finite support. An equivalence class of syndromes is called a superselection sector. A syndrome 
is called topologically trivial iff it is equivalent to the trivial syndrome s(u) = for all u. 

In Section [231 we have constructed closed flexible loops lying in [lll]-bilayers. We shall now need 
to construct such closed loops for all other bilayers. Let t = (t x ,t y ,t z ) be any body-diagonal, see 
Eq. (TjO), and a be any odd integer. Define a plane 

St, Q = k) G A odd : it x +jt y + kt z = a}. 
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For any given syndrome s we define a family of topological Z2 charges as 

s(w) (mod 2). (19) 

Note that the sum above is well-defined since, by definition, s has finite support. 

We will now show that these topological charges 9t, a form a complete set, that is, two syndromes 
s, s' are equivalent, see Definition [31 iff 9t, a (s) — 9t,a{s') for all t and a. By linearity, it suffices to 
prove the following. 

Lemma 3. A syndrome s can be caused by a Pauli operator P with finite support if and only if for 
all odd a and all body- diagonals t, one has 9 a ^{s) = 0. 

To prove the "only if" direction we note that any Pauli operator P with finite support can be 
enclosed by a sufficiently large closed flexible loop lying in any [i]-bilayer. The corresponding string 
operator Wt,a(j) must commute with P since their supports are disjoint. On the other hand, Wt, a {l) 
coincides with the product of generators S v over all plaquettes of the bilayer that are encircled by 
7, see Section I2~3l Hence the commutation between P and Wt, a {l) is controlled by the topological 
charge 9 t<a (s), and we conclude that 9 t , a (s) = 0. To prove the other direction, we manipulate the 
syndrome by operators with finite support to effectively annihilate it. 

Proof. Suppose P G V has finite support and let s be its syndrome. Let us show that 9t >a {s) = 0. 
Since P has finite support, we can take any sufficiently large flexible loop 7 lying in the [t]-bilayer 
£t,a-i U Xt )Q +i such that 7 encloses the support of P. Obviously, the corresponding closed-string 
operator Wt, a {l) commutes with P, since their supports are disjoint. On the other hand, the operator 
W t , a {l) is a product of the stabilizer generators in its interior Int(7), Eq. ( Tl3|) . and thus PWt ta (j) = 
(— l) e *.«( s )Vt / t)Q! (7)P, i. e. the closed-string operator picks up the topological charge of P. Thus the 
topological charge must be zero, 9 t>a (s) = 0. 

Now we prove the other direction, namely if for a syndrome s all charges 9 t a (s) = 0, then s can be 
caused by a Pauli operator P with a finite support. Applying single-qubit X-errors we can shift the 
support of s onto the pair of [010]-planes with y-coordinates j — and j = 1. Applying single-qubit 
F-errors we can shift the support of s in each of these planes onto a pair of adjacent [001]-lines i — 
and i — 1. Since the composition of all X- and F-errors used above has finite support, it does not 
change functions 9 t , a (see the first part of the proof). Hence we can assume that s has support only 
on the union of lines Loo, L01, L10, and Ln, where 

Lab = k) G Kdd ■ i = a,j = b}. 

Suppose that s ^ 0. Let k be the largest integer k such that s(a, b, k) = 1 for some a, b G {0, 1}. Let 
u = (a, b, k). Consider four cases. 

Case 1: (a, b) = (0, 0). Then the plane with t = (—1, —1, 1) contains a single non-zero syndrome 
bit s(u) = 1, u = (0, 0, k). It means that 9 tt k{s) = 1 which is a contradiction. 
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Case 2: (a,b) = (0,1). Then the plane S ti _fc-i with t = (1,-1,-1) contains a single non-zero 

syndrome bit s(u) = 1, u — (0, 1, k). It means that 6 t -k-i{s) = 1 which is a contradiction. 

Case 3: (a,b) = (1,0). Then the plane Y^t-k-i with t = (—1,1,-1) contains a single non-zero 

syndrome but s(u) = 1, u — (1, 0, k). It means that 9 t .-k-i{s) = 1 which is a contradiction. 

Case 4'- («, b) = (1, 1). Then the plane T, t k +2 with t = (1, 1, 1) contains a single non-zero syndrome 

bit s(u) = 1, u = (1, 1, k). It means that Qt,k+2{s) = 1 which is a contradiction. 

Summarizing, we get a contradiction unless s = 0. It means that the original syndrome can be 
caused by a Pauli operator with a finite support. □ 



3.3 Decomposition into dipoles and quadrupoles 

Let s be any syndrome (an excited state) such that the total number of excitations in s is even. 
In this section we show how to decompose s into a combination of dipoles and quadrupoles. To 
illustrate the usefulness of our formalism we show how to create a 'dislocation' in a rigid string that 
can be composed only of quadrupole charges located in the neighborhood of the dislocation. 

Let us first determine the configuration of topological charges, Eq. (1TUI) . corresponding to a single 
quadrupole, a dipole and a monopole, the results are summarized in Fig. [TTJ 
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Figure 11: Given a syndrome s, the corresponding superselection sector can be described by a diagram 
where the black dots indicate non-zero topological Z2 charges Q^ a - Here t labels body-diagonals and 
a is an arbitrary odd integer. Examples of diagrams representing a single quadrupole of type [111], 
a single dipole of type [110], and a single monopole located at site u = (0, —2, 1) are shown. For any 
syndrome s with even (odd) number of excitations, the number of dots on each horizontal line must 
be even (odd). For any syndrome s the number of dots with a = 1 (mod 4) and the number of dots 
with a = — 1 (mod 4) must be even. 

Consider as an example a quadrupole of type t = [111] that consists of four excitations located 
at the sites 

Uq — x, u — y, Wo — z, u + x + y + z (20) 

for some site uq G A even , see Lemma [TJ Let s be the corresponding syndrome. Let a be the odd 
integer such that u 6 S t>a _i. Then 6 ty p{s) = 1 iff (3 = a ±2. Consider now some other body-diagonal 
t' 7^ t. Let a' be the odd integer such that m G £t', a '-i- Then one can easily check that two of the 
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four sites in Eq. (120]) belong to the plane S t / a /_ 2 while the other two sites belong to the plane S 4 / a /. 
Hence 9t',p(s) = for all t' ^ t and all j5. Similar arguments can be made to quadrupoles of all other 
types, see Lemma [TJ We conclude that for any odd integer a and for any body-diagonal t there exists 
a quadrupole which has only 9 t ^_ 2 and 9 tjCe+2 as non-zero topological charges, see Fig. [TTJ 

Similarly we can determine the configuration of topological charges, Eq. (lT9j) . corresponding to 
a single dipole of type h. Consider as an example a dipole of type h = [110] (that can be created 
by a rigid string lying in the [001]-plane) that consists of two excitations located at the sites u — x 
and Uq — y, see Section 12.21 These two sites belong to the same plane £ t)C , for the body-diagonals 
t = [111] and t = [111]. Hence the corresponding topological charges 9t,p are zero for all /3. On the 
other hand, the sites uo — x and u$ — y belong to two different planes £t, Q ±i if t — [111] or t — [111], 
where a is the even integer such that uq G S t a . Hence we have four non-zero topological charges 
9t,a±i where t = [111] or t = [111], see Fig. [Til Note that these are the two body-diagonals which lie 
in the plane orthogonal to h. Similar arguments can be applied to dipoles of all other types. 

Now we are ready to prove that any syndrome s with an even number of excitations can be 
decomposed into dipoles and quadrupoles. First, we note that for any such syndrome 

Em*) = E s ( w ) = °- ( 21 ) 

a u 

Secondly, we have an identity 

E E m*) = E E m*) = o, (22) 

t a=l (mod 4) t a=— 1 (mod 4) 

where t runs over the set of body-diagonals, see Eq. (jl]). Indeed, let u = k) £ A 0( ^, t = [t x t y t z ] 
be some body- diagonal, and u ■ t = it x + jt y + kt z . Since i + j + k is odd, the inner product u ■ t 
is odd for all body-diagonals t. One can easily check that Ylt u ■ t = 0, where the sum runs over all 
body-diagonals t. Hence the number of t such that u-t = 1 (mod 4) must be even. Analogously, the 
number of t such that u ■ t = —1 (mod 4) must be even. It implies Eq. f|22p . 

Let us represent s by a diagram as shown on Fig. [TTJ such that each non-zero charge 9t, a is 
represented by a dot. Combining s with a single quadrupole we can shift any dot along the horizontal 
axis by distance 4, that is, (t, a) — > (t, a ± 4). It allows us to concentrate all dots in the pair of 
columns a = ±1. Since 6 , t _ 1 (s) + 9 tl i(s) = 0, see Eq. f[2Tl) . we are left only with four charges 

9 t (s) = 9 t>1 (s) = 9 t ^ 1 (s) 

labeled by body-diagonals t, see Eq. (TT]). From Eq. (1221 we infer that 

E^ s ) = > 

that is, the number of non-zero charges 9t(s) must be even. Any such syndrome s can be decomposed 
into a single dipole or a pair of dipoles. Indeed, suppose 9 t (s) = 1 for some pair of body-diagonals 
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ti,t2- Then s can be created by a single dipole of type h, where h is the face-diagonal orthogonal 
to t\ and t 2 , see Fig. [11] for an example. In the remaining case, 9t(s) = 1 for all body-diagonals t. 
Then s can be created by a pair of dipoles of types [110] and [110]. Summarizing, any syndrome s 
with even number of excitations can be decomposed (in the non-unique way) into a combination of 
quadrupoles and dipoles. 

3.4 Rigid strings with dislocations 

Let us now consider a pair of parallel rigid strings 71,72 that are shifted with respect to each other 
as shown on Fig. [T2J We would like to connect the two strings together obtaining a rigid string with 
a 'dislocation' and explore what topological charges describe such a dislocation. 

Let W(ji), W("f2) be the string operators corresponding to 71,72, see Section [2T2l Recall that 
W(7i) creates a dipole near each end-point of 71. Assume without loss of generality that the dipole 
located near the dislocation occupies a pair of sites 

= (-1,2,0) and v x = (0,1,0). 

Similarly, W^v) creates a dipole that occupies a pair of sites 

u 2 = (1,0,0) and v 2 = (2,-1,0) 

and another dipole on the opposite end-point of 72. Hence an operator 1^(7172) = W (71) W{j2) 
creates four excitations at sites u\,Vi,U2,V2- Let s be the corresponding syndrome. One can easily 
check that the only non-zero topological charges 6 tya (s) correspond to body-diagonals t = [111], 
t = [111] and a G {—3, — 1, 1, 3}, see the diagram shown on Fig. [131 One can easily check that this 
diagram is equivalent to a pair of quadrupoles of type t = [111] and a pair of quadrupoles of type 
t = [111], see the diagram of a single quadrupole shown on Fig. [TTJ For any t as above one has 
one quadrupole in the [t]-bilayer H t) -2 U £ tiQ and one quadrupole in the [t]-bilayer H t>0 U S t)2 . We 
shall see in Section H] that for a finite lattice with periodic boundary conditions satisfying g = 1 one 
has only one [t]-bilayer for each body-diagonal t (different planes H t , a must be identified due to the 
periodic boundary conditions). Hence for a finite lattice each pair of quadrupoles of type t can be 
paired up with flexible strings and annihilated by the corresponding string operators. Note that if 
the lattice dimensions are pairwise co-prime, a rigid string aligned with a face-diagonal cannot be 
closed directly since the corresponding diagonal fully fills up a two-dimensional plane before it gets 
closed. Therefore creating a dislocation (or several dislocations) might be the only possible way to 
construct a closed rigid string for such lattice geometry. 

4 Encoding of a qubit using closed rigid strings 

Let us remind the reader of some standard notations pertaining to stabilizer codes, see e. g. [T2] . 
Let V = (Xi, Zi, . . . ,X n , Z n ) be the n-qubit Pauli group. Given a Pauli operator P, its weight |P| 
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Figure 12: A pair of parallel shifted rigid strings 71 and 72 of type [110] can be connected into a 
single rigid string 7172 with a dislocation. Double circles indicate locations of excitations associated 
with the dislocation. 
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Figure 13: The diagram shows topological charges associated with a dislocation on a rigid string, 
see Fig. [T2J These charges are equivalent to a pair of quadrupoles of type t = [III] and a pair of 
quadrupoles of type t = [111] located in adjacent [t]-bilayers. 



is defined as the number of qubits on which P acts non-trivially (as X, Y, or Z). A stabilizer code 
can be defined by an abelian subgroup S C V such that —I S. The corresponding codespace C is 
spanned by n qubit states invariant under the action of any element of S. Pauli operators commuting 
with every element of S form the centralizer of S denoted as C(S) = {P G V : PQ = QP VQ G5}. 
The elements of the centralizer preserve the codespace C. For a stabilizer code encoding k logical 
qubits the centralizer can be represented as C(S) = (S, A 1; Z\, . . . , X k , Z k ), where Xi and Z^ are 
logical Pauli operators on the z-th encoded qubit. In this case dim£ = 2 k . The distance of the code 
S is defined as minimum weight of a logical operator, d = mirip & c(S)\s \P\- 

If the code encodes more than one qubit, some of them may be protected better than the others. 
In this case it may be advantageous to use only a subset of k' = k — m 'good' logical qubits to 
encode information. We then consider the gauge group Q = (S,Xi, Z±, . . . ,X m , Z m ) generated by 
stabilizers and the logical operators acting on the remaining m 'bad' logical qubits. Pauli operators 
acting non-trivial on the good logical qubits are elements of C(S)\Q. Hence the distance of the code 
in which quantum information is encoded into the 'good' subsystem is dg = minp 6C ( 5 )\ g \P\. 

In the next subsections, we assume that the lattice has periodic boundary conditions, 

A = Z L x Z L x Z L , 

J_J X Ljy Lj Z J 
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where L a = 2p a for some integers p x ,p y} p z - In addition, we shall focus on the special case when 
p x ,Py,Pz are odd and pairwise co-prime, 

gcd(p x ,p y ) = gcd(p y ,p z ) = gcd(p z ,p x ) = 1, (23) 

p a (mod 2) = 1, a = x,y,z. (24) 

The purpose of the co-primality constraint Eq. (1231) is to take advantage of the rigid strings present 
in the model to maximize the code distance, see Section 14.31 The constraint Eq. (TS4J) is introduced 
to simplify classification of logical operators and is not essential. 

4.1 Half- filled membrane operators 

In this section we construct a complete set of logical operators for the code S. Under conditions 
Eq. ( 1231) . Theorem H] implies that the code S has k = 4 logical qubits. One can choose the logical 
Pauli operators on the four encoded qubits using half-filled membrane operators which are supported 
on four disjoint sublattices of A even . Let a,b,c G {0, 1} with a + b + c (mod 2) = label the four 
qubits. For a qubit labeled by abc, we can find the logical operators 

= n n x aj,k, 

j=b (mod 2) k=c (mod 2) 

= n n y *A*. 

i=a (mod 2) k=c (mod 2) 

= n n z «> (25) 

i=a (mod 2) j=b (mod 2) 

Note that the logical operators a^ bc are supported on a sublattice A abc for which the x,y and z 
coordinates have parity a, b and c respectively. For example, the four logical a x ahc operators lie in 
one of two [100]-planes characterized by fixed x-coordinate % — a — or i — a — 1. For such fixed 
[100]-plane, the choice of b (and thus c) determines the parity of the y-coordinate j and z coordinate 
k, see Fig. HU A simple inspection shows that the operators a" bc commute with all generators S v 
and thus they are elements of the centralizer C(S). 

It is straightforward to see that any pair of logical operators <r" fec and o^, b , c , commute when 
abc 7^ a'b'c' since they have disjoint supports. Consider the commutation between, say, o x abc and a v abc 
and notice that they overlap on p z qubits. Since p z is odd, Eq. ( 12"4"|) . these logical operators thus 
anti-commute on an odd number of qubits, hence they anti-commute. The choice for odd p x and p y 
similarly ensures the proper commutation rules between a x abc and a z abc etc. 

It is important to note that one can translate half-filled membrane operators by multiplication 
with the stabilizer generators. For example, one can shift the logical operator <7g 00 to the parallel 
[100]-plane k) = (2, j = (mod 2), k — (mod 2)) by multiplying it with stabilizer generators at 
the sites k) = (1, (mod 2), (mod 2)). So we can translate any half-filled membrane operator 
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Figure 14: (Color Online) A half-filled membrane operator a z lying in some [001]-plane is the product 
of Pauli Z's over all qubits which have a fixed parity of the x- and ^/-coordinates (the qubits indicated 
by double circles). 



lying in a plane orthogonal to the basis vector a by any vector 2ka for integer k and get an equivalent 
logical operator. 

Although the logical Pauli operators introduced in Eq. (|25|) have a membrane-like geometry, it 
does not mean that any logical operator of the code S has membrane-like geometry as well. Indeed, 
it may be possible to transform the half-filled plane operators or some combination of them into a 
logical operator with a smaller support (e.g. string-like) by multiplying them with the stabilizers. In 
fact, we will see in the next section that some combination of these logical operators are obtainable 
as closed strings winding around the lattice, see Lemma HI 



4.2 Logical operators associated with closed flexible strings 

Recall that flexible strings are confined to [t]-bilayers U S tiQ!+1 , where t is a body-diagonal and 

Ht,a is a [t]-plane that includes all sites u = (i,j,k) satisfying t x i + t y j + t z k = a. Note that for 
periodic boundary conditions and lattice dimensions satisfying Eq. (1231) the sum t x i + t y j + t z k is 
defined only modulo gcd(L x , L y , L z ) = 2gcd(p x , p y ,p z ) — 2g — 2. Hence the entire 3D lattice can be 
considered as a single [t]-bilayer embedded into the 3-torus and folded into a closed surface. One can 
easily check that any such bilayer has topology of the two-dimensional torus^l It follows that flexible 
strings can be closed in a topologically non-trivial way by winding around the torus. We would like 
to characterize the corresponding string operators in terms of the logical operators Eq. (|25|) . Let us 
define the following combination of the logical operators: 

^1 = ^000^011^101^110' ^2 = ^ooo^oi 1^101 ^iio- (26) 

2 For example, one can check that any [t]-bilayer is a closed orientable surface with Euler characteristic x = since 
it admits covering by the regular hexagonal lattice. It follows that this surface has genus g — 1, that is, it has the 
topology of the 2-torus. 
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In this section we will show that Z\, Z 2l and Z\Z 2 are the only logical operators that correspond to 
closed flexible strings winding around the lattice. On the other hand, elements of S correspond to 
contractible flexible loops, see Section 1231 Hence the subgroup of C(S) generated by closed (trivial 
or non-trivial) flexible strings can be described as 

g'=(S,Z u Z 2 ). (27) 

Lemma 4. The group Q' is the group generated by all string operators ^(7) associated with closed 
flexible strings 7. 

Proof. First, recall that stabilizer generators can be regarded as homologically trivial flexible loops, 
see Section 1231 In addition, it is clear that W(j) G C(S) for any closed flexible string 7. Let us first 
prove that 

W(i)a: bc = (-1)^^(7), (28) 

where (w x , w y , w z ) G H^S 1 x S 1 x S 1 , Z 2 ) represents the homological class of 7. In other words, w a is 
the parity of the number of times the a = x,y,z coordinate of the string 7 winds around Zj, a . Using 
the lattice symmetries, it suffices to consider the case when a = z and 7 lies in the [lll]-bilayer. 
Let M abc be the support of a z bc (a half-filled membrane) and n be the [001]-plane containing M abc . 
Let 7 = (uq, Ui, . . . , Ut = uq). Recall that 7 is a sequence of steps X x , X y , X z with alternating signs, 
see the definition of A" in Eq. ([7]) in Section 12.31 Without loss of generality, the string never makes 
two consecutive steps that undo each other such that Uj + 2 = Uj for some j. We shall say that a site 
Uj is a crossing point iff Uj G 7r and Uj+i = Uj ± \ z (and hence both Uj, Uj+i belong to the plane 
7r). Note that exactly one of the sites Uj,Uj+i belongs to the half-filled membrane M a b c - Since the 
steps (uj-\,Uj) and (uj+i,Uj+2) must be ±\ x or ±\ y , the string operator W(j) anti-commutes with 
a z ahc at exactly one of the sites Uj,Uj + i. Taking into account that the sites and Uj + % are on the 
opposite sides of tt, we infer that the number of sites at which Wi^y) anti-commutes with a z bc and 
the winding number w z have the same parity. 

Since each pair of consecutive steps u — > u + A" — preserves the sum of the coordinates i+j + k, 
the winding numbers of any closed flexible string must obey the constraint 

w x + w y + w z = (mod 2). (29) 

If all w a are even, then Eq. (|2"81) implies that W(j) commutes with all logical operators, that is, 
W{^f) G S and hence ^(7) is a trivial loop. When W(j) G C(S)\S for some closed flexible string 
7, exactly two of the numbers w x ,w y ,w z must be odd. Consider for example the case when w x ,w y 
are odd and w z is even. Combining Eqs. (I28II29P we infer that W{^) anti-commutes with all logical 
operators a x bc , a v abc and commutes with all logical operators a z abc . There is only one operator in 
C{S)\S which has this property and this is Z 2 = Ylabc^abc- Similarly, when w x is even and w y ,w z is 
odd we infer that W{^) should be Z\ = Ylabc^abc- Thus the group generated by all closed flexible 
strings W(j) equals the group (S, Z\, Z 2 ). Note that the logical action of non-trivially closed flexible 
strings is only determined by a topological property, namely their winding numbers. □ 
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One may ask how the minimum weight of nontrivially closed flexible strings scales as a function 
of n, where n is the total number of qubits. Let 'j m in be the shortest closed flexible string which 
is topologically non-trivial (i.e. non-contractible to a point). As was mentioned above, one can 
consider the entire 3D lattice A as a single [t]-bilayer £ 4 associated with some body-diagonal t which 
is embedded into the 3-torus and folded into a 2-torus. Note that S t has area of order n since it 
contains n qubits covering the bilayer with density of order 1. Then Loewner's torus inequality implies 
that the shortest non-contractible loop on the torus E t has length \^ m i n \ < c^fn for some constant 
c. It means that the minimum weight of non-trivial closed flexible strings cannot grow faster than 
Cy/n. In fact, one can show that for a particular sequence of lattice dimensions p x ,Py,Pz = 6(723) 
satisfying Eqs. f l2"3"f!H|) we also have a matching lower bound |7 m j„| = Q(\/n), see [T8] . 



4.3 Logical operators associated with closed rigid strings 

Consider as an example a rigid string associated with the face-diagonal h = [110], that is, a sequence 
of sites 7 = (uo, u%, . . . , u m ) where u i+ i = Ui + h for all i, see Definition CD In order to get a closed 
string, u m = uq, one must have m = (mod L x ) and m = (mod L y ). The smallest positive m 
satisfying these conditions is the least common multiple of L x and L y , that is, 

m = hih = Apxp y = 2px p (30) 

gcd(L x ,L y ) 2gcd(p x , p y ) xy 

Thus the string 7 has to wind p x times around the y-axis (p y times around the x-axis) before it gets 
closed. Such a string 7 completely fills up the xy-p\ane that contains the site u , that is, 7 is actually 
a closed two-dimensional membrane. The corresponding string operator W^) can be expressed as 
W{^) = o"ooo^iio or W(t) = ^101^011 depending on whether 7 lies in [001]-plane with even and odd 
^-coordinate. This operator has support on a fully-filled membrane since it involves all qubits lying 
in some [001]-plane. One can similarly construct logical operators associated with other closed rigid 
strings 7. Such operators generate a subgroup 

/"> // i — z — z — z — z — y — y — y — y — x — x — x — x \ / o i \ 

y = Wooo^iiO' ^loi^oii' ^ooo^ioi' ^on^iiO' ^ooo^oii' ^no^ioi/- 

By comparing Eq. ff27|) and Eq. fl3Tl) we conclude that string operators associated with closed rigid 
and closed flexible strings pairwise commute, that is, PQ = QP for all P G Q' and Q G Q" . It allows 
us to choose a subsystem encoding with two logical qubits and two gauge qubits such that string 
operators associated with closed flexible strings act only on the gauge qubits. For example, let us 
complement Z\, Z2 defined in Eq. ( 126|) to a complete set of logical operators as 

^101 ^110 

0"i m 0"i 1 n 



x x 


— ^OOO^Oll^llO 


Zi 


= ^"000^011 


x 2 


= ^101 


z 2 


= ^000^011 


x 3 


= ^000^011 


Z-i 


= ^000^110 


x 4 


= ^OOO^llO 


z. 


= ^OOO^Oll 
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One can check that Xi, Zj obey the commutation rules of the Pauli operators on four qubits. We can 
now regard qubits 1, 2 as gauge qubits (which do not encode any information) while the remaining 
qubits 3,4 are the logical qubits of the code. The corresponding subsystem code has a gauge group 

Q = (S, Zi, Z 2l X 1 ,X 2 ). 

Lemma H] implies that string operators associated with closed flexible strings act only on the gauge 
qubits. Interestingly, the logical operators A3 and Z3 correspond to fully-filled membranes, while the 
logical operators X4 and Z4 correspond to a pair of half-filled membranes occupying two adjacent 
[100]- and [001]-planes respectively. 

In order to find the distance dg of this subsystem code one has to minimize the weight of the logical 
operators X 3 , A 4 , Z 3 , Z 4 by multiplying them with stabilizers and gauge operators. This distance 
must obey the general upper bound dg = 0(L 2 ), where L is the smallest of the lattice dimensions, see 
Theorem 1* in [12]. Computing the distance dg leads to a highly non-trivial optimization problem. To 
illustrate the potential difficulties let us assume that the lattice dimensions obey one extra constraint 
Py — Px + 2. In this case one can construct a closed rigid string 7 lying in some [001]-plane that 
winds around the lattice only once by creating a single dislocation, see Fig. [T21 As was explained 
in Section 13. 4[ the quadrupoles located in the neighborhood of the dislocations can be annihilated 
by pairing them up with flexible strings 5. The total length of these flexible strings is \8\ = 0(\/n), 
see Section 14.21 It allows us to construct a logical operator W = W{ r y)W{S) whose support has 
geometry of a string-net — a collection of rigid and flexible strings interconnected with each other. 
The total length of the string- net 7 U 5 grows at most as 0{^Jn). Note that the logical operator W 
anti-commutes with the fully- filled membrane operator A3. Indeed, we can assume that the fully- 
filled membrane corresponding to A 3 crosses the rigid string 7 at exactly one qubit and hence W 
anti-commutes with A3. On the other hand, A 3 commutes with the flexible string operators W(5), 
see Eq. ( 1281) . assuming that the end-points of 5 are sufficiently far from the membrane corresponding 
to A3. It follows that W anti-commutes with A3 and thus W cannot be in the gauge group, that is, 

dg < \W\ =0(y/n). 

In the general case the minimum-length string-nets corresponding to closed rigid strings are more 
difficult to analyze. Whether the distance obeys the upper bound dg = 0(y/n) for an arbitrary 
choice of the dimensions p x ,P y ,Pz is an interesting open problem. 

5 Discussion and open problems 

We would like to remark that our model suggests natural generalizations of models other than the 
surface code to the third dimension. Consider for example the qudit surface codes [19] defined on 
a lattice with qudits with prime-dimension d which realizes a 2D gauge theory. To get a three- 
dimensional extension, we can take the lattice A = A even U A D dd defined above and for every u G A dd 
we define a generator S u = X u -% X~^_~ Y u -y ^u+ y Z u -z Z~+ & + h.c. where X and Z are generalized 
Pauli-operators with X d = Z d = I and X' 1 = X f and Z a X b = exp(2ma ■ b/d)X b Z a . One can easily 
check that S u are mutually-commuting (use [A ® B, <g> A>] = which holds for any generalized 
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Pauli operators A, B). An open question is also how one can extend the quantum double models 
based on non-Abelian gauge groups [I] to the third dimension in a similar fashion. 

A natural question is whether the studied 3D model exhibits thermal stability such that the 
relaxation time of error corrected logical operators grows as a function of the lattice size. We believe 
that similar to the 2D toric code, topological order in our model is lost at any non-zero temperature. 
Consider for example a simplified noise model with only Z errors. In this case one can consider the 
dynamics of excitations created in any [001]-plane separately. In any given [001]-plane one may expect 
a gas (whose density depend on the temperature) of dipole strings whose ends freely oscillate along 
diagonal lines. It is favorable energetically for dipole strings to combine to former longer oscillating 
strings which diagonally wind over the lattice and eventually can become the logical operator Z% 
(the fully-filled membrane lying in the [001]-plane). In other words, one expects that the thermal 
behavior of some of these excitations is like those of a ID Ising model which has been folded into a 
2D plane. 
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A Ground state degeneracy 

In this Appendix we prove Theorem [U It suffices to show that the stabilizer group S has n — k 
independent generators with k = 4g and S does not contain —I. 

Let us begin by introducing some notations. Let M be any finite set and B = {0, 1}. A set of 
bit strings of length \M\ in which bits are labeled by elements of M will be denoted B(M). Given 
a string s G B(M) and u G M let s u G B be the bit corresponding to u. For any s,t G B(M) let 
sQ)t G B(M) be the bit-wise XOR of s and t. A subset C C B(M) is called a (binary) linear subspace 
iff s © t G C for any s,t G C. The dimension of a linear subspace C defined over the binary field 
will be denoted dimC. Note that \C\ = 2 dimC for any subspace C. A subspace C that includes only 
constant strings s = (00 ... 0) and s = (11 ... 1) will be referred to as a repetition code on M. 

Recall that g = gcd(p x ,p y ,p z ). We shall begin in Section IA~T1 by proving Theorem[T]in the special 
case g — 1 since in this case the proof is much simpler. It will allow us to introduce necessary 
machinery for the general proof, see Sections IA.2fA.3[ 
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A.l The special case g = 1 

Define a linear subspace C = C(p x ,p y ,p z ) C £>(A 0(W ) describing linear dependencies among the 
generators S^, 

C(p JB ,p y ,Pz) = {teB(A« ia ) : fj S*r = ±I}. 

By abuse of notations we shall omit the dependence on p x ,p y , Pz whenever it is not important. Clearly 
the stabilizer group S has n — dim C independent generators and thus the number of logical qubits is 

k = dimC 

Using the explicit form of generators S u we infer that t G C iff t satisfies parity checks 

t u -x © tu+x © tu—y © tu+y = (32) 
tu—x © tu+x © tu—z © tu+z = (33) 
tu—y © tu+y © tu—z © ^m+z = (34) 

for every site u G A even . 

Lemma 5. For any t G C one has 

t u — t u +2gx — t u +2gy = t u+ 2 gz f or M U G h-odd- (35) 

If g = 1 then the converse is true, that is, any string t G B(h dd) satisfying Eq. |23]] belongs to C. 

In the special case g = 1 Lemma implies that the subspace C can be regarded as a product 
of four repetition codes defined on the four non-overlapping sublattices A 10 o, A 010 , A 00 i, and A m , 
where 

A a bc = {(i, j, k) G A : i = a (mod 2), j — b (mod 2), k = c (mod 2)}. 

Since the product of generators S u over each of these sublattices gives +/ we conclude that dim C = 4 
and —I^S thus proving Theorem [TJ In the rest of the section we prove Lemma [51 

Proof. By taking linear combinations of the parity checks Eqs. (I32H34P one easily gets 

tu-hx © t u +hx © t u -hy © t u +hy = (36) 
t u -hx © tu+hx © t u -hz © t u +hz = (37) 
tu-hy © tu+hy © t u —hz © t u +hz = (38) 

for any non-negative integer ft,. Let us try to choose h such that the above parity checks become 
equivalent to the ones in Eq. fl35l) . We shall use the following well-known fact. 
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Proposition 1 (Chinese Remainder Theorem). A system of equations 

h = ai (mod ni) 
h = a 2 (mod n 2 ) 

has a solution h iff a\ = a 2 (mod n), where n = gcd(ni, n 2 ). 

Let g a p = gcd(p a ,pp). Note that gcd(L Q ,,L /3 ) = 2g a/3 . By Proposition [T] we can choose integers 
h, h' such that 

f 2h = (modLj,) ^ f 2/i' = (mod L z ) 

\ 2h = 2g xy (mod L x ) \ 2h! = 2g xz (mod L x ) 

Using these h and h! in resp. the parity checks Eq. fl36l) and Eq. fl37|) we get 

tu tu+2g xz xi 

for all u G A 0( ^. By Bezout's identity, the smallest positive integer obtained as an integer linear 
combination of 2g xy and 2g xz is equal to gcd{2g xy ,2g xz ) = 2gcd(g xy , g xz ) = 2g, that is, we get 
t u = t u+2gx for all u G A 0( ^. The remaining parity checks in Eq. (|35|) are obtained analogously. In 
the special case g — 1 the parity checks Eq. ( J35l) specify a product of four repetition codes defined on 
the non-overlapping sublattices Aioo, A io, A oi, and A m . One can easily check that such repetition 
codes satisfy all parity checks Eqs. (13211341) . Hence in this case Eq. ( 1351) is a necessary and sufficient 
condition for t G C. □ 



A. 2 Computing the number of logical qubits for g > 1 

We can use Lemma [5]to establish an isomorphism between C(p x ,p y ,p z ) and C(g,g,g) as follows. Let 
A' = Z 2g x Z, 2g x Z, 2g . Given a string t G B(A odd ) obeying Eq. ( 135]) . let t' G B(A' odd ) be a string defined 
through t' a a = tij t k, where a = i (mod 2g), (3 = j (mod 2g), and 7 = k (mod 2g). Lemma [5] tells 
us that t G C(p x ,p y ,p z ) iff Eq. ( )35|) holds and t' G C(g,g,g). Hence in order to prove the identity 
k = 4g it suffices to show that 

dimC(g,g,g) = 4g. (39) 

Let us begin by considering a two-dimensional version of this problem. Consider a 2D square lattice 
Q = Z 2g x Z, 2g . A site u = (i,j) G Q is called even (odd) iff i + j is even (odd). Let Q e ven and fl dd 
be the even and odd sublattices. Define a linear subspace 

C(g, g) = {te B(Q odd ) ■. t l+1J © u_ ltj © t iJ+1 © = for all (ij) e n even }. (40) 

We claim that 

timC{g,g) = 2g. (41) 
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Indeed, considering t = {tij} as a binary matrix of size 2g x 2g (with zero entries at all even cells), 
one can put an arbitrary binary strings at the rows i = 0,1 since there are no parity checks that have 
support only in these rows. Using the parity checks Eq. (|4*0I) one can uniquely fill up the remaining 
rows i = 2, 3, . . . , 2g — 1. Furthermore, using the generalized parity checks Eq. (l3T)j) with h = g one 
arrives at t i+2g j = tij, that is, expressing the row % = in terms of the rows % = 2g — 2,2g — 1 
one will always come back to the original binary string {to.j}- Similarly, expressing the row % = 1 in 
terms of the rows i — 2g — 1, one will always come back to the original binary string {tij}. Hence 
t G C(g,g) is uniquely determined by 2g bits t j, (0, j) G fi dd, and (1, j) G These bits can 

assume arbitrary values. It proves Eq. ( )4T|) . The above arguments also show that we can ignore the 
parity checks centered at the first and the last rows, that is, C(g,g) can be specified by the parity 
checks 

t i+ i,j = t iJ+1 © tij-! © tj-ij, j G Z 2g , 1 < i < 2g - 2, G Vt even . (42) 

Remark: For the later use we show examples of strings t G C(g,g) on Fig. [151 One can easily check 
that these strings and their translations generate the entire subspace C(g,g). 
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Figure 15: Examples of a string t G C(g,g) for g = 3,4 (only non-zero bits are shown). Note that 
these strings have a single non-zero bit in the first pair of rows. It follows that the entire subspace 
C(g,g) is spanned by horizontal and vertical translations of the shown string. 

Now we can easily prove Eq. f )39p . Consider any string t G C(g,g,g). We note that for any site 
u G A even the parity checks Eqs. (I32|33|34p are linearly dependent, namely, 

(t u -x © tu+x © tu-y © tu+y) © (tu-x © tu+x © tu-z © tu+z) © {t-u—y © tu+y © tu-z © tu+z) = 0. 

Using the parity checks Eq. f )32|) we conclude that the restriction of t onto any xy-plane must belong 
to C(g,g). Hence any xy-plane potentially contributes 2g to the dimension of C(g,g,g). We claim 
that t is actually uniquely determined by the restriction of t onto the pair of x?/-planes (i, j, k = 0) 
and (i,j,k = 1). Indeed, applying the parity checks t u+ % = t u _z © t u _ x © t u+x , or the equivalent 
parity checks t u+ ^ = t u _^ © t u _ y © t u+yi one can uniquely extend t to the remaining xy-planes 
k — 2, 3, . . . , 2g — 1. Moreover, this extension automatically satisfies the parity checks Eq. (|32|) in 
every xy-plane k = 2, 3, . . . , 2g — 1. Indeed, consider any pair of sites v = u + z. Suppose u belongs 
to the /c-th x?/-plane and we have already checked that the parity checks Eq. ( 132]) are satisfied in all 
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x?/-planes 0,1, ... ,k. Let us check that 

tv—x © tv+x © t'u-y © ^u+y = 0- (43) 

Indeed, applying Eq. (134")) to the sites u ± x we get 
Similarly, applying Eq. (|33|) to the sites ti ± y we get 

tv±y ^n+a;±j/ © t u —x±y © tu—z±y- (45) 

Adding together Eqs. ( I44"f4"5]) we arrive at 

ty—x © ^f+x © t v —y © t v -\-y ^u— z— £ © ^tt— £+:r © ^u—z—y © ^tt— z+y 

since we assumed that the parity checks Eq. (I32I) are satisfied in the (k — l)-th xy-plane. 

To summarize, we have shown that C(g,g,g) is isomorphic to the direct sum of two copies of 
C(g,g) associated with the pair of xy-planes (i,j,k = 0) and (i,j,k = 1). Hence Eq. (|39|) follows 
from Eq. (T4T]) . 

A. 3 Proving that — / ^ S for the case g > 1 

It remains to prove that —I ^ 5. Consider any string i 6 C such that 

S(t)= J] St=e(t)I, e(t) = ±l. (46) 

«eA odd 

Since the generators S u pairwise commute and = I we conclude that e(t) is a homomorphism 
from C to Z 2 , that is, e(t © £') = e(t)e(t') for any t, t' G C. Thus it suffices to check that e(t) = +1 
only for basis vectors t G C. To this end we define a two-dimensional version of the generators in 
Eq. (TjQ), namely, 

where u G A odd . Note that S' u is the restriction of the 3D generator S u onto the xy- plane that contains 
the site u. Define also an operator 

S'U) := J] (S' u y*. (47) 

ueA odd 

For any xy-plane Q a generator S u with u ^ Q acts on qubits of only by / or Z. Hence the 
cancelation in Eq. fj4"6"]) is possible only if S'(t) is a Pauli operator of Z-type. More specifically, one 
must have 

S'(t) = e(t) J] Z? 
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for some string s G B(A even ). Since each Z v can be only obtained from the product X V Y V or Y V X V , 
we conclude that 

To determine the phase factor e(t) let us choose a specific ordering of generators in Eq. (I47j) such 
that generators S' u with odd x-coordinate are always on the left of generators with even x-coordinate. 
(Since the generators S' u pairwise commute, the order does not matter.) Then any qubit v with even 
x-coordinate can only produce Z v through a multiplication X V Y V . Any qubit v with odd x-coordinate 
can only produce Z v through a multiplication Y V X V . Hence we arrive at 

where n even and n odd is the number of sites v G A even with s v — 1 such that v has even and odd 
x-coordinate respectively. If we sum Eq. ( 148)) over all sites v in A even with odd x-coordinate, we see 
that n odd must be even. Hence we arrive at 

e(t) = i n —n+n odd = ^ ^ 

where \s\ is the Hamming weight of s. Using the identity \a ffi b\ = a + b — 2a6 which holds for any 
a, 6 G {0, 1} and Eq. ( l4"8"j) one can easily get 

uGAnoUAioi ueAoooUAon 

^ ^ 4|t n | 2 ^ ^ tvfu+2x t u tu+2y- 
mGAooiUAoio neAooiUAoio 

Now it follows from Eq. ( )49l) that 

e(t) = (-l) /{t) , /(*) = Yl tutu + 2 * + (mod 2). 

"SAooiUAoio 

Since e is a homomorphism, it suffices to check that f(t) = for all basis vectors of C. We shall use 
the isomorphism between C(p x ,p y ,p z ) and C(g,g,g) described in the beginning of Section IA.2I Let 
t' G C(g,g,g) be the image of t under this isomorphism. Using Eq. ( 1351) one can easily get 

fit) = P -^f(t>). 

g 6 

Here fit') is calculated for generators S u defined on the lattice A' = 7j 2g x ^2g x ^2g- Hence it suffices 
to check that /(£') = for all basis vectors t' G C(g,g,g). As was mentioned in Section IA.21 the 
restriction of C(g,g,g) onto any xy-plane Q = 7L 2g x 7L 2g coincides with the subspace C(g,g), see 
Eq. ( l40l) . Let r G C(g,g) be the basis vector defined on Fig. [T5l Since r is symmetric under 90° 
rotations of the lattice, one can easily get 

^ r u r u+2x + r u r u+2y = (mod 2). 

u=(2i,j)eQ odd 

Hence any xy-plane yields even contribution to fit'), that is, fit') is even. Thus f(t) is even for all 
t G C(p x ,p y ,p z ) which implies e(t) = 1. We have proved that —I ^ S for g > 1. 
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B Proof of zero-temperature stability 



In this Appendix we explicitly prove that the Hamiltonian Eq. (J2j) defines a stable phase of matter 
at zero temperature, that is, the spectral gap of the Hamiltonian does not close in the presence 
of weak local perturbations and its ground state degeneracy is not lifted up to exponentially small 
corrections. Specifically, we will check the sufficient conditions for stability derived in [17]. We also 
show that the model Eq. (j2J) has a unique ground state if defined on an infinite lattice. 

The following lemma proves two statements about Pauli operators P whose support can be 
bounded by a sufficiently small box B. First, if such an operator P commutes with all generators S u 
then P must be in the stabilizer group (up to a phase factor). Secondly, if P is a stabilizer then one 
can express P as a product of generators S u using only those generators supported inside the box B. 
The first statement implies that the stabilizer code S = (S u , u G A D ^) has a macroscopic distance 
(growing at least linearly with the smallest of the lattice dimensions). It provides the first stability 
condition called TQO-1 in Ref. [T7J. The second statement is equivalent to the second stability 
condition called TQO-2 in Ref. [T7J, see Lemma 2.1 in the above reference. As was shown in [T7j . 
conditions TQO-1, 2 together are sufficient for zero-temperature stability. 

Lemma 6. Suppose a Pauli operator P commutes with all generators S u , u G A C( ^ and suppose the 
support of P can be bounded by a box B of size l x x l y x l z with l a < L a — 3 for all a — x,y, z. Then 
P is a stabilizer up to a phase factor. Moreover, there exists a subset of generators M C B D A Q dd 
such that P = YiueM ^ u U P ^° a phase factor. 

Proof. Let I = m&x(l x ,l y ,l z ). We prove the lemma using induction in I. The base of induction is 
/ = 2. Note that a box of size 2x2x2 contains four even sites. Using the translational symmetry 
we can assume that P acts non-trivially on at most four qubits located at sites u\ = (0,0,0), 
u 2 = (1, 1, 0), Us — (1, 0, 1), and u 4 = (0, 1, 1). For any site Wj one can choose a triple of generators 
S v acting on Ui by X, Y, and Z, and acting trivially on the remaining three sites Uj, j ^ i. Hence 
the commutativity PS V = S V P implies that P = I up to a phase factor. 

Let us now prove the step of induction. Assume without loss of generality that l z > l x ,l y . We 
can also assume that l z > 3 since otherwise we have the base of induction. We shall construct a 
stabilizer S G S such that PS has support in a box of size l x x l y x (l z — 1). This stabilizer S will 
only use generators whose support is contained in B. Indeed, let F be the upper face of B. The 
commutation PS V = S V P implies that P u G {/, Z} for all u G F. In addition, P u = I iff u lies on an 
edge of F. Define the stabilizer S as 

S= J] S u . z . (50) 

u£F:P u =Z 

Note that S has support in B and PS acts trivially on F. Hence PS has support in a box of size 
l x x l y x (l z — 1). It proves the step of induction. □ 

A simple corollary of the lemma is that the model Eq. (J2]) has unique ground state if defined on 
an infinite lattice, A = Z x Z x Z. Indeed, define a stabilizer group S* generated by finite products 
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of generators S u , u G A m- For any finite subset M C A even let S(M) be the subgroup including all 
elements S G S* whose support is contained in M. We note that —I ^ S(M) since otherwise one 
would have —I G S for a finite lattice of sufficiently large size which contradicts to results obtained 
in Appendix [A] Hence S(M) can be regarded as a stabilizer code by itself. Define a mixed state pu 
proportional to the projector onto the codespace of S(M), 

PeS(M) 

Note that for any pair of subsets M' C M the inclusion S(M r ) C S(M) implies pm' = T^m\m' Pm- 
Hence any pair of states pm, Pk with overlapping supports have consistent marginal states on the 
intersection M(~)K. A collection of such states p = {pm} associated with all finite subsets M C A even 
defines a quantum state of the entire lattice. By construction one has S u pM = Pm for any generator 
S u whose support is contained in M. Hence p minimizes every term in the Hamiltonian Eq. ([21) 
thus being the ground state of the model. We claim that this state is pure, namely, p cannot be 
represented as a mixture of two different quantum states. Indeed, suppose there exists a pair of 
quantum states r = {r M } and r] = {t]m} such that 

p M = (1/2) (tm + Vm) (51) 

for all finite M C A even . Let us show that this is possible only if p M = r M = r\ M . Indeed, since pu is 
proportional to the projector onto the 5(M)-invariant subspace, we conclude that the range of tm 
is spanned by (S(M)-invariant states, that is, 

Pt m = t m P = t m for any P G S(M). (52) 

Choose any P G S* and let K = Supp(P) U M. Consider the identity Pt^P = tk which follows 
from Eq. ([52]) . Taking a partial trace over K\M we get Pm^mPm = t m, where Pm is the projection 
of P onto M. Thus commutes with projections of stabilizer onto M. Hence any Pauli operator 
P in the Pauli expansion of tm commutes with projections of stabilizer onto M and thus P G «S(M) 
up to a phase factor, see Lemma [6j We arrive at 

PeS(M) 

for some real coefficients a(P). Using Eq. ( )52|) we conclude that a(P) = a(I) for all P G <S(M), that 
is, tm = Pm- The same argument shows that t]m = Pm- Thus we proved that r = rj = p, i.e., p is a 
pure state. 
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